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ABSTRACT 


We  describe  the  Starfire  Optical  Range  Atmospheric  Monitor  (SAM).  SAM  is  designed  to  measure  the 
Fried  coherence  length  r0,  the  Greenwood  frequency  fc-  and  the  isoplanatic  angle  even  in  relatively  poor 
conditions  (ro  ~  1  cm)  with  high  accuracy  ~  10%  with  a  single  frame  of  data.  In  addition,  since  SAM  is 
based  on  a  Shack-Hartmann  wave  front  sensor,  more  sophisticated  statistics  about,  the  wave  front  and  the 
atmosphere  should  be  able  to  be  recovered  from  the  data  with  additional  processing. 


1.  INTRODUCTION 


Many  different  methods  of  measuring  atmospheric  turbulence  have  been  used  over  the  years  ([1],  [2],  (3), 
i4)).  The  Starfire  Optical  Range  Atmospheric  Monitor  (SAM)  was  designed  to  help  us  understand  the 
characteristics  of  the  atmosphere  when  turbulence  conditions  are  poor,  such  jus  during  the  clay,  and  to  also 
gather  “higher  order”  information  for  wavefront  statistics  not  normally  calculated  by  typical  r0  meters. 

In  order  to  collect  even  under  severe  conditions,  SAM  incorporates  a  Shack-Hartmann  wave  front  sensor 
with  small  (2  cm)  subpertures.  This  enables  SAM  to  measure  values  of  ?•(,  somewhat  smaller  than  2  cm. 
Moreover,  we  should  be  able  to  estimate  even  smaller  values  when  we  enable  the  operation  of  the  full  4x4 
array  of  camera  pixels  associated  with  each  subaperture. 

The  collection  of  “higher-order”  information  is  enabled  by  our  Shack-Hartmann  wave  front  sensor  (WFS). 
Because  of  the  large  number  of  subapertures,  we  collect  a  statistically  significant  sample  of  tilts  during  each 
frame,  enabling  rapid  estimation  of  r o,  and  hence  the  ability  to  detect  rapid  changes  in  ro-  In  addition,  using 
a  Shack-Hartmann  WFS  enables  us  to  measure  not  just,  a  few  separations  as  are  sampled  by  a  differential 
image  motion  sensor,  but  many  separations  simultaneously.  We  intend  to  leverage  this  information  to  make 
measurements  of  deviations  from  the  Kolmogorov  model,  and  also  to  enable  the  measurement  of  different 
atmospheric  layers  as  described  below. 


2.  THE  INSTRUMENT 


The  opt  ical  tube  assembly  (OTA)  for  SAM  includes  a  stock  lf>”  Ritchey-Chretien  telescope  from  RC  Optical 
Systems  (RCOS)  .  The  RCOS  telescope  has  a  secondary  which  obscures  the  inner  31%  of  the  subapertures. 
The  OTA  is  mounted  on  a  Paramount  ME  mount.  The  pointing  accuracy  for  the  Paramount  ME  ranges 
from  10  to  30  arcseconds  after  performing  a  mount,  model. 

SAM  is  located  in  a  small  10.5  foot  dome  1G  meters  above  the  top  of  the  hill  that  is  referred  to  as  “Mount 
Fugate”.  SAM  is  situated  a  few  meters  from  the  1.5  in  and  3.5  m  site  telescopes.  Fig.  1  is  a  picture  of  SAM. 


Fig.  1.  Dome  location  (left),  SAM  (right). 


3.  WAVE  FRONT  SENSOR  OPTICS 


Liglil  from  the  RCOS  telescope  is  passed  to  a  custom  optics  assembly  mounted  at  the  Cassegrain  back  plate 
of  the  telescope  (a  schematic  of  the  optics  is  shown  in  Fig.  2).  First,  a  square  field  stop  at  the  focus  of  the 
RCOS  eliminates  light  which  might  fall  onto  adjacent  subapertures  (see  Fig.  3).  The  field  stop  is  slightly 
undersized,  in  that  it  obscures  approximately  one  pixel  around  each  sulmperture.  The  outer  ring  of  pixels 
is  used  as  a  "guard  band",  and  is  ignored  in  all  calculations.  This  limits  the  magnitude  of  the  tilt  that  the 
instrument  can  accurately  measure,  and  hence  the  minimum  value  of  r0  that  it  can  measure. 

After  passing  through  the  field  stop,  the  light  is  collimated  and  passes  through  a  filter  wheel  holding  six 
1-inch  round  filters.  500,  (iOQ,  650  and  700  tun  long  pass  filters  are  available;  the  best  Idler  is  selected  bast'd 
on  the  star  selection  to  maximize  the  SNR  without  saturating  the  detector. 

Next  an  image  of  the  pupil  is  formed  on  t he  lenslet  array,  which  forms  20  x  20  subapertures  across  t he  pupil 
of  the  telescope.  Another  lens  relays  the  images  formed  by  the  lenslet  array  to  the  CCD;  each  lenslet  image 
is  formed  on  4  x  4  pixel  elements  on  the  CCD.  The  held  of  view  for  one  pixel  is  40/irad,  while  the  size  of 
each  subaporture  is  .02m. 

The  CCD  itself  is  mounted  to  tin1  optics  wedge.  The  camera  and  associated  electronics  are  an  off  the  shelf 
configuration  from  SciMeasure.  The  CCD  is  an  80  x  80  pixel  CCD-30  from  E2V  with  a  read  noise  of  9 
electrons  and  a  peak  quantum  efficiency  of  approximately  59%.  The  camera  system  is  capable  of  doing  full 
frame  reads  at  a  rate  of  2kHz;  in  older  to  maint  ain  a  sufficiently  high  signal  to  noise  we  typically  operate  at 
1kHz,  capturing  1000  frames  in  one  second.  These  frames  are  saved  to  HDF5  files  for  analysis.  Typically,  a 
mean  value  of  r o  is  calculated  using  2000  frames  every  45  seconds  when  SAM  is  being  operated  as  a  facility 
instrument . 
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Fig.  2.  Schematic  of  the  Optics.  Light,  travels  from  the  top  to  the  bottom. 


4.  ESTIMATING  /■„ 


Fig.  3.  Depiction  of  CCD  pixels  for  two  subapertures. 


For  the  first  incaniat  ion  of  SAM.  we  estimate  the  till  using  the  differential  tilt  method,  hut  instead  of  using  a 
pair  or  a  small  number  of  subapert tires,  we  use  many  pairs  of  subapert  ares  across  the  whole  Shack-llart  maim 
wave  front  sensor.  For  now  we  estimate  the  value  of  rt (  from  the  mean  differential  lilt  as  discussed  in  [5], 
and  we  only  calculate  the  value  of  i\,  based  on  adjacent  subapertures.  We  differentiate  between  4  different 
estimates  in  this  paper.  First,  we  take  adjacent  apertures  in  either  the  x  direction  or  the  y  direction.  Second, 
we  can  take  the  difference  in  the  gradient  either  parallel  to,  or  perpendicular  to,  the  vector  separating  the 
two  subaperturcs  in  question.  As  we  shall  see  in  the  results  section,  these  4  estimates  all  agree  with  each 


other,  as  they  ought  duo  to  thu  theory  presented  in  [5],  which  assumes  Kolmogorov  phase  statistics.  We  first 
calculate  the  gradients  from  the  energy  imbalances 


Or 


nh+ h - h ~ h 

1 1  +  I.i  -f  I‘2  +  I  \ 


n  was  calibrated  by  injecting  known  tilts  in  the  lab.  Given  these  angles,  we  calculate  the  differential  tilts 


~  °n,yf)  and  rxjj  A  n  =  ((0Aa  -  61)<X)'2), 

where  (•)  represents  averaging  over  several  subapertures.  Finally  these  quantities  are  converted  into  ro.i 
and  t'o.n  using 

o\  =  0.1119(//r0)6/3(A//)2  and  a\ f  =  0.0484(//r0)5/3(A//)2, 
the  formulae  found  in  [5]. 

The  agreement  of  these  1  statistics  (see  the  next  section)  represents  a  very  simple  confirmation  of  the  Kol¬ 
mogorov  theory.  In  future  work,  we  will  be  using  other  statistics  to  further  test  the  Kolmogorov  hypothesis. 
For  example,  by  looking  at  subapertures  beyond  the  adjacent  ones,  we  can  test  the  5/3  power  law  of  the 
structure  function,  and  perhaps  detect  the  inner  scale  in  some  conditions. 


5.  VALIDATION 


We  have  validated  both  the  algorithms  and  basic  design  of  SAM  against  simulations,  and  validated  the  basic 
functionality  of  SAM  against  another  ro  instrument. 
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Fig.  4.  Simulation  of  SAM  measurements 


Ill  order  to  validate  the  optical  and  algorithmic  design  of  SAM,  we  created  a  wave  optics  simulation  of 
SAM  in  the  ACS  wave  optics  simulation  system.  We  propagated  through  a  variety  of  atmospheres,  collected 
simulated  data,  and  ran  the  data  through  our  algorithms.  The  results  are  shown  in  Fig.  4.  The  blue  circles 
represent  atmospheres  with  It)  Hz  Greenwood  frequencies,  the  red  squares  200  Hz  Greenwood  frequencies. 
We  conclude  that  the  software  and  hardware  systems,  as  designed,  should  measure  i'(>  down  to  small  values 
of  /■(>  even  when  fa  is  quite  large. 

We  also  compared  SAM  ro  measurements  from  another  ro  meter  deployed  nearby,  the  "mini-SAM"  instru¬ 
ment.  mini-SAM  estimates  ro  from  the  MTF  of  the  image  formed  during  a  short  exposure.  mini-SAM  is 
known  to  have  biases  for  large  values  of  r o.  In  Fig.  J»  we  compare  observations  with  the  two  instruments:  we 
see  good  agreement. 


Fig.  5.  Comparison  of  SAM  to  mini-SAM. 


(i.  ESTIMATING  0o 


The  isoplanatic  angle,  as  dclined  by  Fried  [(>].  weighs  the  optical  turbulence  C; ;  with  a  c  ,/:1  altitude  depen¬ 
dence: 


Walters  and  Bradford  [7]  estimate  0q  from  measurements  of  the  normalized  variance  of  intensity 
aperture: 


0,,=  Bcoss/r,(0) 


-3/5 


across  an 


B  refers  to  a  scaling  function  determined  by  the  size  of  the  apert  ure  calculated  in  Walters  [7].  a  is  the  signal 
variance  of  the  irradiance  across  the  aperture,  while  S  is  the  average  signal.  is  calculated  by  summing 
the  signal  of  each  frame,  calculating  the  standard  deviation  across  2000  frames,  and  dividing  the  standard 
deviation  by  the  average  signal  for  those  2000  frames. 


SAMs  larger  apert lire  (.  Ini)  follows  a  quadratic  z~  dependence.  Fig.  <i.  To  achieve  l lie  '  scaling  l lie  SAM 
aperture  was  virtually  apodized  by  .!).  creating  a  smaller  aperture. 

SAM  0{)  measurements  with  the  apozided  filter  were  once  again  compared  to  mini-SAM,  Fig.  7.  We  see 
good  agreement  between  the  instruments. 


Optical  Weight  of  Tuibulenca  vs  Attitude 


Fig.  (i.  C,y  weights  with  altitude  for  SAM  and  Do 
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F’ig.  7.  ft()  lucasurenients,  SAM  and  mini-SAM. 


7.  FROZEN  FLOW  AND  TOMOGRAPHY 


Shock  iiiid  Spilliir  [.'{]  demonsl ruled  t hat  the  st at istics  of  individual  at niosplierie  layers  could  lie  teasi'd  out  of 
wave  front  sensor  data.  They  used  data  from  the  1.5m  telescope  at  SOIL  and  emphasized  the  validation  of 
the  frozen  How  hypothesis.  Since  SAM  is  also  based  on  a  wave  front  sensor,  it  should  be  possible  to  repeat 
and  extend  this  work. 

The  crux  of  that  work  is  the  observation  that  if  we  regard  the  measured  pupil  phases  (or  gradients)  as  a 
three  dimensional  c  ube  of  data  <!>(:r.y.t),  and  also  assume  that  the  frozen  How  hypothesis  is  valid,  then 

d>l(.r  -  t'jT.  i/  -  v,jT.  t  t)  =  0/(.r,  i/.  y) 

for  each  frozen  layer  /  for  all  times  r  so  long  as  r  is  short  compared  to  tin*  time  at  for  which  the  layer  remains 
frozen.  If  we  take  the  three  dimensional  autocorrelation  of  the  data  cube,  this  relation  implies  that  there  is 
a  strong  autocorrelation  along  a  trajectory 

(A  r,  Ay,  At)  =  (uxt,vvt,t). 

These  trajectories  are  visible  as  "streaks"  in  the  autocorrelation.  As  the  frozen  component  breaks  down, 
the  intensity  of  the  peak  decays  at  larger  t.  If  we  add  more  layers,  it  is  easy  to  see  that  each  contributes 
its  own  streak  following  its  own  trajectory,  under  the  assumption  that  each  of  the  random  phase  screens  is 
independent  of  i  he  other  phase  screens.  Finally,  although  we  have  discussed  only  the  phase  screens  0(.r,  y,  y), 
other  quantities  which  co-move  with  the  phase  screens,  including  gradient,  and  intensity  fluctuations,  also 
follow  the  same  general  behavior. 
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Fig.  8.  Intensity  auto  correlation  of  data. 


Using  early  SAM  data,  we  have  reproduced  the  analysis  not  with  gradients  as  in  [Sj.  but  with  the  intensity 
fluctuations,  to  produce  the  data  displayed  in  Fig.  8.  We  see  that  the  autocorrelation  behaves  as  expected: 
evidently  the  wind  was  blowing  in  the  positive  x  and  positive y  din'd. ion.  (The  signal  is  rather  weak  partially 
because  of  the  rudimentary  signal  processing  employed,  we  believe  that  further  signal  processing  will  tease 
t  la*  signal  out .) 

Future  work  will  attempt  to  use  a  combination  of  gradient  and  intensity  data  to  estimate  distances  to  the 
layers,  I  he  speeds  and  directions  of  the  layers,  and  the  relat  ive  intensil  ies  of  the  layers  on  a  second  bv  second 
basis. 


8.  SAMPLE  DATA 


In  Fig.  9  \v<*  display  some  sample  data  obtained  at  SOU  with  SAM.  At  the  top  of  the  figure  we  display  tin* 
value  of  t'(|.  sealed  to  zenith  and  fiOOutu.  for  a  day.  using  both  r  and  u  separ.it ions  and  ±  and  estimators. 
Note  the  excellent  agreement  between  these  measures.  In  the  bottom  left  we  plot  the  standard  deviation  of 
the  measured  values  of  r o-  The  standard  deviation  is  calculated  as  the  root  mean  square  of  the  to  values 
which  occur  during  a  2000  frame  capture,  with  a  separate  r o  calculation  being  performed  on  each  frame. 
Note  that  the  error  in  the  estimate  in  n>  from  a  sint/lr  frame  of  data  is  on  the  order  of  10%.  Finally,  in 
the  lower  left,  we  display  the  signal  to  noise  ratio.  Note  that  the  signal  to  noise  ratio  is  quite  high.  The 
jumps  in  the  SNR  correspond  to  changes  in  stars,  while  the  more  gradual  changes,  like  that  after  l.'i  GMT, 
correspond  to  the  sun  rising. 


s,d(ro)  SNR 


Fig.  9.  Sample  observations. 


9.  CONCLUSIONS 


In  this  paper  we  have  described  tin*  Starhre  Optical  Range  Atmospheric  Monitor.  SAM  is  currently  able  to 
estimate  /•„  accurately  ±10%  with  a  single  frame  of  data,  and  estimate  ()u  from  scintillation  values.  SAM  can 
operate  during  both  day  and  night,  and  operate  during  relatively  severe  (r o  %  Farm)  conditions.  Finally, 
we  are  beginning  to  use  the  Shark-Hurt  maun  wave  front  sensor  data  to  tease  apart  the  layers  of  turbulence 
in  the  atmosphere,  and  estimate  higher  order  statistics  of  the  wavefront. 
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